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We examine the stability of the nonextremal D1-D5-P black hole solutions. In particular, we look 
for the appearance of a superradiant instability for the spinning black holes but we find no evidence 
of such an instability. We compare this situation with that for the smooth soliton geometries, which 
» , were recently observed to suffer from an ergoregion instability, and consider the implications for the 

■ fuzzball proposal. 
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I. INTRODUCTION 

Despite being well-understood classically, black holes still pose a number of unanswered questions at the quantum 
level, such as the information paradox. A radical new approach to describe stringy black holes, now known as the 
'fuzzbair proposal, was first suggested some years ago by Mathur and collaborators They advocated that the 
microstates underlying a black hole are individually described by horizon-free geometries and that the black hole 
geometry only emerges in a coarse-grained description which 'averages' over the e SBH microstate geometries. In this 
approach, the effective horizon of a black hole appears as a surface at a radius where the individual microstate 
geometries start to 'differ appreciably' from one another. Therefore quantum gravity effects are not confined close to 
the black hole singularity but rather the entire interior of the black hole is 'filled' by fluctuating geometries - hence 
£C) • the nomenclature: the 'fuzzball' description of black holes. 

Finding evidence to support this conjecture is not easy though, especially because it entails finding families of 
horizon-free geometries sufficiently extensive to describe all of the e SBH microstates for a given black hole. Most of the 
studies to date focus on supersymmetric configurations, namely the BPS D1-D5 system the BPS D1-D5-P system 
0,13, and the BPS D1-D5-P-KK system [|. However, if the fuzzball proposal is to be useful, it must also be extended 
to non-supersymmetric systems. This then poses the extremely difficult problem of finding non-supersymmetric 
(smooth) horizon-free geometries that can be associated with the microstates of a non-BPS black hole. So far, the 
only known solutions in this class are those of Jejjala, Madden, Ross and Titchener [||, hereafter referred to as JMaRT 
solitons. The JMaRT solutions comprise a five-parameter family of D1-D5-P non-supersymmetric smooth geometries 
which are asymptotically flat. These solutions may be parameterized by the Dl-brane and D5-brane charges, the 
(asymptotic) radius of the internal circle with Kaluza-Klein momentum, and by two integers m and n which fix the 
remaining physical parameters. These integers also determine a spectral flow in the CFT which allows the underlying 
microstate to be identified. For m = n+ 1, the JMaRT solitons reduce to supersymmetric solutions found previously 
in @, 0] . The geometry of these solutions was also recently examined in more detail in 0] • 

In a previous paper [8| , we have shown that the non-supersymmetric JMaRT solitons are classically unstable against 
an ergoregion instability @. This kind of instability is generic to geometries with an ergoregion but no horizon, as was 
first noticed in [9] . Thus this instability should be a robust feature of any smooth horizon- free geometry corresponding 
to a non-BPS microstate with angular momentum. At first sight then, this seems to pose a challenge for the fuzzball 
proposal: if horizon-free non-supersymmetric solitons are expected to describe the microstates of non-BPS black 
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holes, then the nonextremal D1-D5-P black hole which results from averaging over an ensemble of such solutions 
should be unstable against an analogous instability. The purpose of the present work is to address this issue. While 
the presence of an event horizon eliminates the possibility of an ergoregion instab ility , there is, however, an obvious 
candidate instability in the black hole case: the superradiant instability [ToL [Til. Il2j . 

In general, spinning nonextremal black holes will exhibit superradiant scattering, whereby an incident wave packet 
can be reflected with a stronger amplitude. Superradiance by itself does not provide a classical instability, but an 
instability can arise if the waves are reflected back and forth. Stated in other words, a superradiant instability is 
present if there are bound states subjected to superradiance. Some examples of black hole systems unstable against 
such mechanism are: i) the black hole bomb where an artificial mirror surrounds a Kerr black hole 11311 ; ii) a massive 
scalar field in a Kerr background [T3| ; hi) small Kerr-AdS black holes [l5[ ; iv) rotating black strings [lOUlll [l2| . This 
instability seems to be the natural extension, to black holes, of the ergoregion instability found in [8|. Therefore, 
according to the fuzzball proposal, one might expect that the black hole family of D1-D5-P non-supersymmetric 
geometries should be superradiantly unstable - however, we find they are not! Nevertheless, as we will discuss, this 
does not present a sharp contradiction with the fuzzball proposal. 

The remainder of this paper is organized as follows: In Section [TTJ we briefly review some of the relevant properties 
of the D1-D5-P family of supergravity solutions. In Sections Mil and IIV1 we study minimally coupled scalar waves in 
the general D1-D5-P background. We write the wave equation in a Schrodinger form and perform an extensive search 
for unstable free scalar modes in the D1-D5-P black hole geometry. We find no unstable modes. The same conclusion 
applies for a massive scalar field. In Section |Vl we discuss our results and consider some of their implications for 
the fuzzball proposal. In particular, we discuss possible ways the fuzzball proposal can consistently incorporate our 
results. The paper contains two appendices: In Appendix^ we show that the wave equation for the general D1-D5-P 
system written in Section IIIII reduces to the form presented in [H, Q when we restrict our analysis to the JMaRT 
solitons. In Appendix [Bl we give the explicit expression for the Schrodinger potentials introduced in section [TV] 



II. PROPERTIES OF THE D1-D5-P FAMILY OF SOLUTIONS 



The black hole and solitonic configurations considered here are solutions of the type lib supergravity equations. The 
solutions all carry the three charges of the D1-D5-P system and so are expected to give a strong-coupling description 
of different (ensembles of) microstates of this system. The supergravity solution is comprised of a (ten-dimensional) 
metric and also a nontrivial dilaton and RR two-form potential. The system is compactified to five dimensions 
on M 4 x S 1 with the D5-branes wrapping the full internal space and the Dl-branes and KK-momentum on the 
distinguished S 1 . The other component of the compactification M 4 is a Ricci-flat fouf- manifold, which we take to 
be either a four-torus T 4 or K2>. The notation is best understood by considering the construction of these solutions. 
One begins with the general solutions of [1, [H, [l8| which contain eight parameters: a mass parameter, M; spin 
parameters in two orthogonal planes, ai,a 2 (which we assume are non- negative without loss of generality); three 
boost parameters, 5i, £5, which fix the Dl-brane, 1 D5-brane and KK-momentum charges, respectively; the radius 
of the S 1 , R; the volume of the four-manifold, V4 (which plays no role in the following). The geometry is described 
by the six-dimensional line element which is given below (see also Equation (2.12) of |y|) and which is parameterized 
by a time coordinate t; a radial coordinate r; three angular coordinates 9, <fi, tp; and the coordinate on the S , y. 

The (ten-dimensional) dilaton and two-form RR gauge potential which support the D1-D5-P configuration are [H, [f| 

e 2 * 10 =H 1 /H B , (1) 

_ f 2i Mcos 2 6 . , . , 

C y ' = = [(a2Cis 5 Cp - aiSic 5 Sp)dt + [aiSic 5 c p - a 2 c 1 s^s p )dy\ A dip 

Hi 



M sin 2 



H 



[(aids 5 c p - a 2 siC5S p )dt + (a2SiCsc p - aicis^s p )dy] A > 



-M^£L dt Ady — (r 2 + a 2 + Ms 2 ) cos 2 9dip A dcj>. (2) 

Hi Hi 

However, these will only play an ancillary role in the present discussion. Central to our analysis will be the six- 
dimensional geometry consisting of the noncompact space, as well as the y-circle. The contravariant components of 



1 In the following, only the asymptotic Dl-brane charge Qi appearing in the supergravity fields will be relevant. Of course, with M 4 = K3, 
Ql tx Ni — iVg where Ni and iVs are the numbers of constituent Dl- and D5-branes, respectively, comprising the system [Tsj . Ref. \2(j \ 
describes how this technical point produces an interesting physical effect for the D1-D5-P black holes considered here. 
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the string-frame metric are 2 

g u = c 2 J l + 28,0,9® + s 2 p 9™ , g ty = s p c p {g rt + s «) + (c 2 + S 2 )/^ , g™ = s 2 /' + 2s pCp g® + c 2 5 ^ 
g ** = C J* + Spff 5* , g t * = c p g i *+8 p g* l >, g»* = sj* + c p g^, 9^ = 3^+0^, 



jrr 



1 9 9 1 „^ 1 / 1 (r 2 + a 2 )(a 2 -a 2 )-Ma 2 



1 Ma lQ2 ^ = 1 / 1 + (r 2 + a 2 ) (a 2 - t 



1 Maio 2 ^_ 1 / 1 (r 2 + a;;)(a 2 .-a 2 )- Mq| \ 



where 3 



.9" = ~7i^= ( /(r) + M + Af.s 2 + Ms 2 + 

ty _ _ 1 M 2 sis 5 cic 5 aiQ 2 t -0 _ 1 Mcic 5 a 2 (r 2 + a 2 ) ^ _ 1 Afcic 5 ai(r 2 + of) 

5 V#i#5 #( r ) ' 9 VHvH~5 Sir) '-' y%H~ 5 9{r) 

g» = ( f(r) + Ms\ + Msl + M2s ^ 2 + a2 + a2 - M) 



y/llJh V 9(r) 

y s 1 Ms lSbai {r 2 +a 2 - M) w 1 M Sl s 5 a 2 (r 2 + a 2 - M) 
.9 = / - - t~\ ' 9 = — / ~ - n • ( 4 ) 



We are using the notation a = coshJi and Sj = sinh<5i. Throughout these expressions, we also use the functions: 

g{r) = (r 2 + a\)(r 2 + a 2 ) - Mr 2 = (r 2 - r 2 + ){r 2 - r 2 ) , (5) 
f(r) = r 2 + a\ sin 2 8 + aj cos 2 9, Hi(r) = f(r) + Ms 2 , with i = 1, 5 . (6) 

Without loss of generality, we will assume a\ > a 2 > in the following. 

Depending on the values of the parameters, this solution can describe a black hole or a naked curvature singularity. 
However, it was also realized that in a third parameter regime [6|, this geometry corresponds to a smooth soliton - 
denoted the JMaRT soliton in [8( - or a conical singularity. The best way to identify each one of these branches of 
solutions is to look at the g rr component of the general metric, which is proportional to the function g(r), given in 
eq. ((5]) above. The roots of g(r), r + and r_, are given by @ 



r± = \ (M ~a\- of) ± \ yj (M a 2 a 2 ) 2 - 4a 2 a 2 , (7) 

and they are real whenever \M — a\ — a|| > 2aia 2 . We can naturally divide the general family of solutions into three 
branches, namely Q: 

M < (oi - a 2 f ^>r 2 + <0, (JMaRT branch) , 

(ai - a 2 ) 2 < M < (oi + a 2 ) 2 =>■ r 2 ^ K , (Naked singularity branch) , (8) 
M > (oi + a 2 ) 2 > , (Black hole branch) . 

In the first case the system can describe either a smooth soliton or a conical singularity. (In Appendix [XJ we present 
the constraints which must be satisfied to produce a smooth geometry.) In the second case, the function g(r) does 
not have real roots and the system describes a naked singularity, with curvature singularities where Hi(r) vanishes. 
Finally, in the third case, the system describes a black hole with outer horizon at r 2 = r 2 ^ and inner horizon at r 2 = r 2 . 
The upper bound of the first branch, M 2 = {a\ — a 2 ) 2 (for which r 2 = r 2 = — aia 2 ), includes as special cases the full 
set of supersymmetric smooth solitons @, @|. The lower bound of the third branch, M — [a\ + a 2 ) 2 , corresponds to 
an extremal black hole with r 2 , = r 2 = aia 2 . The supersymmetric limit of the above three-charge system corresponds 



2 The boosted coordinates (i, y) are related to the unboosted coordinates (t, y) by t = tc p — ys p , y = yc p — ts p - see Appendix A of [f|. 

3 Eqs. (A. 2), (A.3) and (A.4) of H have a typo and must be multiplied by minus one. 
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to taking the limit M — > and 6i — > oo, while keeping the other parameters fixed, including the conserved charges 
Qi = MsiCj. One gets, in the non-singular case, the black hole solutions of [2l[ or the supersymmetric solitons of 
0,1 -seel. 

In this paper, the key feature of interest in this geometry is the ergosphere. In particular, we will be considering 
instabilities that arise due to the existence of an ergoregion. To verify the presence of the ergoregion, one may take 
the norm of the Killing vector V = dt yielding 

f - M c 2 

9abV V = _ . (9) 

This result shows that V — dt becomes space-like for f(r) < MCp and thus one would conclude that an ergosphere 

appears at f(r) = Ale 2 . However, in the supersymmetric limit, this norm @ becomes | | 2 = — (/ — Q p )/ y H1H5. 
So using this measure, one arrives at the counter-intuitive conclusion that the ergoregion persists even for the su- 
persymmetric backgrounds. The resolution of this puzzle is evident in the discussion of section 6.2 of [||. The key 
point is that the present geometry 'rotates' along the internal y-direction, as well as along the angles (j> and if). Hence 
for the purposes of defining the ergosphere, we have a continuous family of asymptotically time-like Killing vectors: 
V = d t + v y d y with \v y \ < l. 4 Now one can push the position of the ergosphere to smaller radii by adjusting the free 
parameter v y . In particular, the position seems to be minimized by the choice v y — tanh 8 P , for which V = df (i.e., 
we have undone the boost of footnote ^ and 

9 abV a V b = -^L. (10) 

Further we find that this choice of V matches the Killing vector arising from the square of the covariantly constant 
Killing spinor appearing in the BPS limit. Of course, the ergoregion now disappears in this supersymmetric limit where 
M — > and / — > r 2 . A perspective relevant for the following discussion is that in the supersymmetric backgrounds, V 
provides a globally time-like Killing vector which ensures that there is a rotating or 'boosted' frame where all energies 
can be defined to be positive. 



III. SEPARATION OF THE WAVE EQUATION IN THE GENERAL D1-D5-P SYSTEM 

The JMaRT solitons were found to be classically unstable in ref. Q. What we wish to do now is apply a similar 
stability analysis to three-charged black hole geometry to ascertain whether or not these solutions are also unstable. 
A stability analysis starts by disturbing slightly the system and then letting it evolve freely. Stable systems return 
to their original 'position' whereas in unstable ones, the perturbations will grow without bound in time. Ideally one 
would like to consider perturbations by any of the type lib supergravity fields, e.g., metric, dilaton or RR two-form 
perturbations. In Sections |HT] and HVj we will begin by considering only free scalar perturbations that obey the 
Klein-Gordon equation. These perturbations are expected to capture the essential physics while at the same time 
simplifying enormously the calculations. However, as we discuss in Section |V1 our results for the free scalar apply 
directly to certain supergravity fields. 

So we first consider the Klein-Gordon equation for a massless scalar field propagating in a general three-charge 
geometry. The results of this section are valid for the full range of the parameters defining the configuration space. 
Let us note that our strategy will be to solve the wave equation in the six-dimensional background formed by the 
non-compact directions and the distinguished S 1 circle. That is, we limit our analysis so that the directions along 
internal manifold M 4 play no role. This amounts to having a wave with no excitations along these internal directions. 
We must mention that a related analysis of this equation in the effective five-dimensional geometry coming from also 
reducing on the S 1 circle first appeared in [l7|. 5 

We must precede our analysis of the wave equation with the observation that for the D1-D5-P backgrounds, the 
effective six-dimensional dilaton vanishes: The ten-dimensional dilaton is given in eq. ([T]) while the string- frame metric 



4 We note that \v y \ = 1 would yield an asymptotically null Killing vector. Of course, the norm of the angular Killing vectors 8^ and cL, 
diverges asymptotically and so we cannot consider further linear combinations by including either of these vectors. 

5 The present analysis is made in six, rather than five, dimensions because the Kaluza-Klein momentum along the S 1 plays an important 
role, as discussed at the end of that section IIVI Further the coordinates in [l7j | are adapted to calculating the absorption cross-section 
of the black hole rather than producing the Schrodinger form of section IIVI 
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on the internal space M 4 takes the form 



H 



Gab = [ w J Gab = e * 1Q Gab ' 



(11) 



where G a b is a fiducial metric with unit four-volume. Hence upon reducing the supergravity action to six dimensions, 
we find the effective six-dimensional dilaton becomes: 



-2*6 — „-2*io 



VdetG^ 



(12) 



It follows then that in the effective six-dimensional theory, there is no distinction between the string-frame and 
Einstein-frame metrics. In particular, the six-dimensional metric has precisely the same components as those of the 
ten-dimensional string- frame metric, given in eq. ©. 

Hence we consider the Klein-Gordon equation in the effective six-dimensional theory 6 coming from the reduction 
of the general three-charge geometry: 



1 d_ 

dx A 



— AB d 



* =0, 



(13) 



where g AB is the six-dimensional metric explicitly written above, in section [TTJ Further g is the determinant of this 
six-dimensional metric, which is given by 



V— — rsindcosdy HiH 5 . 
This is the same equation studied in the stability analysis of Q. Introducing the separation ansatz 



^ = cxp 



— iw-ji- — iX^- + im^ip + im<p4> 
R R 



X (0)h(r) 



and the separation constant A, the wave equation separates. The angular equation is 

,2 „ 2 



1 



sin 26* rf6» 
and the radial equation is 

1 d 

r dr 



sm29 de 



A 



*</> 



cos- 



m 4> , - X 2 2 ■ 2 a , 2 2 n\ 

(a-L sin + a 2 cos 0) 



sm 



R 2 



x = o, 



~g(r) d h 


-Ah + 


r dr 





(^ 2 -A 2 ) 
R 2 



(r 2 + Ms\ + Msl) + {ujc p + Xs p ) 



M 
R 2 



1 



M 2 
M 2 



2 2 2 2 o ,22 2/ 2,2,2 ji,r\l 

c^c^cj ~ 2s 1 s 5 s p cic 5 c p aia2 + s 1 s 5 s p (r + a 1 +a 2 — M)\ 



2uiX- 



R 2 L ~ s p c p L~ c i c 5 r + s i s 5( r + a 1 +a 2 - M)\ + (c- + s p )s 1 s 5 cic 5 a 1 a 2 \ 

+ [cic 5 c p a 2 {r 2 + a\) - s 1 s 5 s p a 1 (r 2 + af - M)l 

R 



ZCOTn^h r / 2 2\ / 2 2 t\1 

— - — M \ c\C^c p a\\r + a 2 ) - SiS 5 Spa 2 {r +a 2 -M)\ 
R 



M 



2 2 2 



a2 -fiT [ _c l u 5°p 
2Xm r , 



O O O / O Q O \ 

2sis 5 s p cic 5 c p aia 2 + s 1 s 5 c p (r + a 1 + a 2 — M) 



(14) 



(15) 



(16) 



(17) 



H tt^ m [cic 5 s p a 2 (r 2 + a\) - s 1 s 5 c p a 1 (r 2 + a\ - M)l 

H 

H — M|cic 5 s p ai(r + a 2 ) — sis 5 c p a 2 (r + a 2 — Af ) I 

R 

- ml [(r 2 + a\)(a\ - a 2 2 ) - Ma 2 ] - m 2 , [{r 2 + a%){al - a\) - Maj] + 2m^m^Maia 2 \h = 



This wave equation is equivalent to that which would be produced in reducing the equation of motion of a massless scalar field which 
is minimally coupled to the Einstein-frame metric (and no other fields) in ten dimensions. 
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This radial wave equation is valid for the the general three-charge geometries, including the black hole solutions and 
the JMaRT solitons. However, as shown in Appendix |A1 when the JMaRT constraints (| Aip ~ (|A2|) are imposed, the 
wave equation can be rewritten in a considerably simplified way, namely as in eq. (|A9|) . 

The angular equation (|16p . plus the appropriate regularity requirements, defines a Sturm-Liouville problem, and the 
solutions are known as higher dimensional spin- weighted spheroidal harmonics [221 ] . We can label the corresponding 
eigenvalues A with an index I, A(u>) = Ai m . m ,(u>) and therefore the wavefunctions form a complete set over the 
integer I. In the general case, the problem at hand consists of two coupled second-order differential equations: given 
some boundary conditions, one has to compute simultaneously both values of u> and A that satisfy these boundary 
conditions. However, for vanishing af we get the (five-dimensional) flat space result, A = 1(1 + 2), and the associated 
angular functions are simply given by Jacobi polynomials [22]. For non-zero, but small u ^ Q% we have 

K = l{l + 2)+oU l ^L^ . (18) 

The integer I is constrained to be I > \m^\ + \m$\. We will always assume af - ^T 2 A <C max(m^,TO^) (with i = 1,2) 
and thus A ~ 1(1 + 2). Making this assumption implies we may neglect the terms proportional to a% in the angular 
equation. But given the way A and u> appear in the radial equation, the corrections to A may not be negligible when 
we determine w. To ensure that setting A — 1(1 + 2) is consistent in both the angular and radial equations, we must 
additionally require [see first line of (fT7|) ]: af <§; max (|r+ + M(s\ + s§)|, Mc^j. 



IV. WAVE EQUATION IN THE SCHRODINGER FORM 

In Q , we have shown that the non-supersymmetric JMaRT solitons are unstable against the ergoregion instability. 
The ingredients for this instability are the existence of an ergoregion in a geometry without horizon. If there is a 
counterpart of the ergoregion instability on the black hole side of the configuration space, it seems likely to be the 
so-called superradiant instability. The ingredients for this instability are the existence of an ergoregion around a 
horizon and the presence of 'bound' states within the superradiant regime [ToL fTTL IT2I H3L [Lll. [TEl] . 

The simplest way to find if the non-supersymmetric black hole solutions of the D1-D5-P system are superradiantly 
unstable is to rewrite the radial wave equation in the form of an effective Schrodinger equation, and to study the 
corresponding Schrodinger potentials. We now apply this approach in the following. 

In [8| we found that the introduction of dimensionless coordinate 

X — Q 9~ ) (19) 



along with a new wavefunction H 



h(x)= 7 =L== H(x). (20) 
transformed the radial wave equation (fTT|) to 

d l H + m l ^ +x) 2 ^ - U +)^4- -U-)H = 0, (21) 

where V, U- and U+ are presented in Q. However, we find that this form is no longer appropriate to study the black 
hole sector, since in this case both U- and U+ take on complex values, sufficiently close to the black hole horizon r + . 
Instead, we introduce the 'tortoise' coordinate 7 

dr g(r) 



7 Note that this coordinate is not appropriate to analyze the M < (ai — ai ) 2 branch of solutions since r 2 ^ < in this case, and thus r* 
becomes complex near the origin r 2 = ri. 
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where g(r) was defined in eq. ([5]) and the new wavefunction $ is given by 

h = r- 3/2 $. (23) 

Then (|17p can be written as a Schrodinger equation, 

d 2 <$> , , 

^-F* = 0, (24) 



with 



(r^(r) - Irff'(r) - ^«/(r)) . (25) 



9 -^d r h 



Here the prime denotes a derivative with respect to r, and W^(r) is defined by writing (fl7|) in the form i<9 r 
W(r) /i = 0. As discussed earlier, we assume Aw 1(1 + 2) independent of uj and therefore V is a quadratic function of 



V = -j(uj-V+)(uj-V-). (26) 

The explicit forms of 7 and of the Schrodinger potentials V± is given in Appendix [B] It can be checked that the 
asymptotic behavior of the potentials is 

lira V± = uj sup , (27) 

r— *r+ 

lim V± = ±|A| . (28) 

r — >oo 

where we defined the superradiant factor, 

w sup = m^Sl^R + m^Sl^R — XSl y . (29) 

Here, Sl$, fi^,, and Sly are, respectively, the angular velocities along 4>, ip and the velocity along the S* 1 given by [23j 

„ a 2jr 2 + _ r\cic b s p + a 1 a 2 s 1 s 5 c p 

(r| + J (r^cic 5 c p + aia 2 SiS5Spj r + cic 5 c p + a 1 a 2 s 1 s 5 s p 

Typical forms of the potentials V± are displayed in Fig. [T] In these plots, the 'allowed' regions where the solutions 
have an oscillatory behavior are those where uj is above V+ or below V-, i.e., above or below both curves V±. In those 
intervals where uj is in between the curves of V+ and V— (forbidden regions), the solutions have a real exponential 
behavior. From these plots we will infer the stability of the system. We impose only ingoing- wave boundary conditions 
near the horizon and outgoing at infinity. This means that the asymptotic behavior of the solutions of ((23|) -([24 | . is 

( ) [ ( r 2 _ r 2 = e -feln(r 2 -r=) ; as r ^ r+ , 

where we have defined 

M (ricic 5 Cp + aia 2 sis 5 s p ) , . ,„„. 

ro = — 7-2 51 w _ u ^p> ■ ( 32 ) 

2Rr+(ri_ — r_) 

When the frequency of the wave is such that w is negative, 

\u\ < \uj sup \ , (33) 

one is in the superradiant regime, and the amplitude of any incident wave is amplified upon scattering from the black 
hole. 

The simplest way to appreciate this result is through the following argument (24|: 8 From (fTS)) and (|3ip one has that 
at the horizon the wave solution behaves as \&(t, r)\ ~ e —*ut/R e -f07la(r -r + ) _ -jv^g ph ase velocity of the wave is 



See section II. A of [3] for an alternative argument. 




FIG. 1: Typical Schrodinger potentials for the black hole branch, M > (ai + a^) 2 or r\ > 0, and for modes with no KK 
momentum, A = 0. Figure a) refers to the case > 0, while figure b) corresponds to case < 0. Superradiant modes 
(dashed line) with \u>\ < \ui aup \ do exist but there are no trapped bound states and thus the geometry is not afflicted by the 
superradiant instability. 



then Dph oc — — . Now, the value of this phase velocity can be positive or negative depending on the value of u (when 
we fix the other parameters), so one might question if the first line of (|31|) really describes always an ingoing wave. 
What is relevant for the discussion is not the phase velocity, but the group velocity of the waves. The normalized 
group velocity, u gr , at the horizon is t> gr = 2(r\ - rly+a- 1 = -l. 9 This is a negative value that signals that 

the near-horizon wave solution in (|31[) always represents an ingoing wave independently of the value of w, and thus we 
have the correct physical boundary condition. However, note that in the superradiant regime l|33p . the phase velocity 
is positive and so waves appear as outgoing to an inertial observer at spatial infinity. Therefore energy is in fact being 
extracted from the black hole. 




FIG. 2: Examples of potentials that would correspond to black hole geometries afflicted by the superradiant instability. Unstable 
modes are superradiant modes with \u>\ < \uj sup \ that are also bound states (dashed line). Case a) corresponds to A = while 
case 6) corresponds to A 7^ 0. 

The relevance of superradiance in the present context is that it may give rise to instabilities. The only extra 
ingredient one needs for that to happen is to somehow 'trap' the waves near the horizon. An artificial way of doing 



9 Note that for M > (ai + (L2) 2 one has > 0,10.2- Moreover, cf = 1 + s 2 > sf . 
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TABLE I: Some examples of modes that have the Schrodinger potentials plotted in Fig. [3] The black hole geometry is 
described by the parameters (ai = 32 , a2 = 16, M = 1.01(ai + w^) 1 ,c\ = 5,cb = 1.517, c p — 5,R = 1). The modes have 
(I — m^p = 10 , 77i0 = 0) and A indicated in the first row of the table. Similar Schrodinger potentials are obtained when m$ is 
switched on. 

this would be to enclose the black hole inside a mirror [l3|. Any initial perturbation will get successively amplified 
near the black hole event horizon and reflected back at the mirror, thus creating an instability. This instability is 
caused by the mirror, which is an artificial wall, but one can also devise natural mirrors. For example, consider a 
massive scalar field [14]. Imagine a wavepacket of the massive field in a distant circular orbit. The gravitational force 
binds the field and keeps it from escaping or radiating away to infinity. But at the event horizon some of the field goes 
down the black hole, and if the frequency of the field is in the superradiant region then the field is amplified. Hence 
the field is amplified at the event horizon while being confined away from infinity. Yet another way to understand this, 
is to think in terms of wave propagation in an effective potential. If the effective potential has a well, then waves get 
'trapped' in the well and amplified by superradiance, thereby triggering an instability. In the case of massive fields on 
a four-dimensional Kerr background, the effective potential indeed has such a well. Consequently, the massive field 
grows exponentially and is unstable. It is the presence of a bound state that simulates the mirror, and so without a 
bound state we should never get an instability. 

It was found in [l(| HH, E3] that introducing KK momentum for a massless field in black string (or brane) geometries 
can be equivalent to having a mass and so can also trigger superradiant instabilities. Hence we pay particular attention 
to KK momentum in our analysis below. 

The strategy here is to look for bound states in the effective potential to ascertain whether or not the geometry is 
superradiantly unstable. We start by considering the case in which the waves have no KK momentum, A = 0. In this 
case, at infinity the potentials go to zero, as indicated by (f2"5| . In Fig.[T]we sketch two typical examples of potentials 
that can occur when A = 0. The specific parameters that yield these plots are black holes with (ax = 32 , a-i = 
16, M = 1.01(ai + a 2 ) 2 , Ci — 5 , c 5 = 1.517, Cp = 5,ii = 1), and modes with (I = 10 , m$ = , A = 0) and = 10 
[Fig. [TJa] and = —10 [Fig. Q]b]. The geometry has superradiant modes that satisfy (|3"3"]) . and examples of these 
are represented by a dashed line in these plots. In Fig.[T]a, an incident wave with frequency tu sup < u> < is reflected 
in the potential V- and returns back to infinity with an amplified amplitude. Similarly, in Fig. [TJb, superradiant 
modes are those with < to < io sup that are reflected in the potential V+. In both cases, the pattern speed of the 
superradiant modes, E^, = uj/(Rm^), is negative. This means that to be amplified, the waves must satisfy (|33[) and, 
in addition, rotate in the same sense as the black hole rotation (the angular velocity of the black hole fl^ is negative - 
see (|30p ). The second necessary condition for the superradiant instability - the existence of bound states - is however 
absent in these plots. We have done an extensive search (by varying the parameters ai , M , Cj , R and I , , m^,) for 
superradiant bound states and have found none with A = 0. For the sake of clarity, a typical example of Schrodinger 
potentials that would allow superradiant bound states is sketched in Fig. [5] a. 

We now turn our attention to modes with KK momentum, i.e., with A ^= 0. The KK momentum provides a potential 
barrier of weight A at infinity [see (|2"8)) ] and thus are the most promising when looking for instabilities. However, not 
even in this case does there seem to be an instability, as shown in Fig. [3) This figure represents the several shapes of 
potentials V± that we can get when A ^ 0. None of these plots has bound states and so it seems that the D1-D5-P 
black holes are not afflicted by the superradiant instability. In Fig. [2]b, we plot an example of a case for which the 
system with A ^ would have superradiant instabilities. As we just stated, we found no such case. In Table Q] we 
give specific exam ples of parameters that have the Schrodinger potentials presented in Fig. [3J 

As discussed in |ll| the absence of bound states in the potential seems to be closely related to the absence of stable 
circular orbits in the geometries, which is a generic feature of higher dimensional spaces [l(J [HI, [Hj|. It is not so 
surprising then to find that higher dimensional rotating black holes are stable against the superradiant mechanism. 

We end this section by considering perturbations of a massive free scalar. Massive scalar fields, with mass [i, obey 
the equation 

-4=d A (^g AB d B V) = ^ , (34) 

with the determinant given by Eq. (fT4|) . It is easy to see that, under the same separation ansatz (|15|) . the wave 
equation will be separable if and only if \f—g itself separates. This happens in particular if Hi = H^, i.e., if 
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FIG. 3: Black hole branch, M > (ai + a^) 2 or r\ > 0, for modes with KK momentum, A / 0. Although the KK momentum 
provides a barrier at infinity, it is not able to create bound states where superradiant unstable modes could eventually live. 

Qi = Qs- 10 I n this case, under the assumption of small mass //, the angular equation is still given by (|16p with the 
angular eigenvalue A defined in (| 18[) . and the radial wave equation (JTTJ) receives now a source term on its right-hand 
side given by fi 2 (r 2 + Ms\)h. We followed the same procedure to turn this equation into a Schrodinger-like ODE. 
Again, we did not find evidence of bound states. Therefore the inclusion of a scalar mass does not seem sufficient to 
make these geometries unstable through the superradiant mechanism. 

V. DISCUSSION 

In this paper, we studied a general family of supergravity solutions for the D1-D5-P system which contains two 
special branches: one with smooth horizon-free geometries (JMaRT solitons) and a black hole branch. In general, 
when rotating, these solutions have an ergoregion. This is absent only in the supersymmetric case (in both branches). 
However, we note that if these backgrounds are regarded as solutions of type I supergravity, some of the ergo-free 
backgrounds are non-supersymmetric (on both branches), e.g., certain solutions are extreme but non-supersymmetric 
black holes [23l. |25||. With an ergoregion, the background exhibits superradiant phenomena, however, superradiant 
scattering per se is harmless. It will only extract a small amount of rotational energy from the geometry and transfer 
it to a wave that transports it to infinity. There are, however, two situations where the addition of extra ingredients 
leads to a catastrophic instability. One is when we have a spacetime geometry with an ergoregion but no horizon. 
This situation generically leads to the so-called ergoregion instability Q in which there are modes which are outgoing 



This condition is also satisfied to a good approximation, if ai/r+ and a^/r^ are very small. 
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at infinity, regular at the origin, and growing unboundedly with time. The negative energy that is stored in the 
crgoregion core, by energy conservation, will then also grow negative without bound. In a previous article, we found 
that the D1-D5-P smooth horizon-free geometries are afflicted by this instability 

The other scenario which produces a catastrophic instability is more well-known. A rotating black hole with an 
crgosphere and with a 'reflecting wall' can lead to multiple superradiant scattering/reflection that extracts rotational 
energy from the black hole without bound. This reflecting wall can be provided by an artificial mirror [HI, by the mass 
of the propagatin g fie ld [l4| , by an asymptotic anti-de Sitter geometry [TEj or by the KK momentum of a black string 
geometry [ToL [TH The key feature, necessary for the activation of this instability, is that the effective potential 
that describes the field propagation in the given background must have a well where bound states can be trapped. The 
D1-D5-P black hole could potentially be afflicted by this superradiant instability. Indeed, it has both an ergoregion 
and KK momentum along the distinguished S 1 circle, which might create the mentioned reflecting boundary. However, 
in our extensive search over the parameters of the solution, we have found that the other ingredient - the potential 
well - is absent. Hence the D1-D5-P black holes do not appear to suffer from any superradiant instability. 

In sections HTT1 and llVl we verified the absence of a supperradiant instability for a minimally coupled scalar field (both 
massless and massive) propagating in the D1-D5-P black hole background. Of course, it would be most interesting to 
verify that the black hole remains stable when we perturb it by the fields of the type lib supergravity theory. Hence, 
we now address the question of to what extent our analysis applies to the supergravity fields. Of course, the same 
discussion applies to our instability analysis of the JMaRT solitons @ which was also explicitly carried out with a 
massless minimally coupled field. 

As discussed at the beginning of section IIII1 we are considering the propagation of scalars in the six-dimensional 
space comprised of the non-compact directions and the distinguished S 1 circle. Hence let us consider the scalars 
appearing in the reduced six-dimensional supergravity. After compactifying the type lib supergravity down to six 
dimensions, the scalars parameterize the moduli space 

Mo = SO(5, n; R)/SO(5) x SO(n) (35) 

up to additional global identifications, where n = 5 or 21 for M 4 = T 4 or K3, respectively (26j. Hence the six- 
dimensional theory contains a total of 5n independent scalars. However, when the D1-D5 string is introduced in the 
six dimensions, several of these scalars interact with the RR two-form ([2} sourced by the string. These 'fixed' scalars 
acquire an effective mass in the AdS^ x S* 3 core of the black hole. Setting aside these fixed scalars, the residual moduli 
space is 

M = 50(4, n; R)/SO{A) x SO{n) (36) 

again up to additional global identifications [261 ] . Hence we are left with An 'minimal' scalars in the six-dimensional 
D1-D5-P black string background. 

In the six-dimensional effective action, the kinetic term for the scalars can be written as [13]: 

(fx^gg^Li^M^L^M 11 , (37) 

where the scalar fields are represented by the (5 + n) x (5 + n) matrix M y taking values in the coset (|35|) . These then 
satisfy the following identities: 
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(38) 

where l d are d x d unit matrices and the superscript T indicates matrix transposition. As discussed in section UTT1 
the six-dimensional dilaton vanishes and hence in eq. (]37p . g^ v corresponds to both the string- frame or Einstein- frame 
metric. The full scalar action also includes couplings to the gauge fields and form fields in six dimensions [271 ] . as well 
as to internal fluxes 28, 29]. As discussed above, in the present background, these interactions are irrelevant for the 
An minimal scalars and so eq. (|37[) compromises the entire effective action for these fields. Now this action ([3"T|) has a 
deceptively simple form and so one might conclude that all of the minimal scalars satisfy the massless Klein-Gordon 
equation (fTB")) . However, in fact, the action is implicitly nonlinear since the scalars take values on the coset (13"6"|) and 
so in general, this conclusion is mistaken. For example, with the torus compactification, i.e., M 4 = T 4 , the equation 
of motion for scalars originating from the internal components of the ten-dimensional metric is 

V 2 G afc - G ac G bd \/ 2 G cd = 0. (39) 
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when the other fields are set to zero. Carefully examining an explicit representation of M [27|, one finds that relatively 
few of the scalars actually satisfy eq. (fTS)) in general. Recall that we must focus on the minimal scalars of and 
then there are just three scalars comprising the antisymmetric tensor parameterizing the 4x4 block just above the 
diagonal in the top-right corner of M [27| . i.e., in the same position as 1 5 appears in L. In the torus compactification, 

(2) 

these three scalars correspond to the internal components of the RR two-form, C ab and beginning with the d = 10 
supergravity action, one can easily demonstrate that these field obey the ordinary massless Klein-Gordon equation in 
six dimensions. The general discussion here confirms that in fact these three fields are the only scalars satisfying this 
simple equation of motion, in general. Of course, there is one case which deserves special attention: Qi = Q5. In this 
case, the background scalars, i.e., the ten-dimensional dilaton ([TJ) and the internal volume (fTTj) . are constant. This 
trivializes the nonlinearities in the scalar field equations and any linear perturbations of the minimal scalars about the 
background satisfy the massless Klein-Gordon equation (fT5)) . For example, for the T 4 compactification, the internal 
moduli may be written as G a b — S a b + h a b with traceless perturbations h a b and the eq. (|39j) reduces to V 2 h a b = 0. 

The lesson we derive from the above discussion is that there are precisely three fields in the six-dimensional 
supergravity whose fluctuations are described by the Klein-Gordon equation (fl~3|) analyzed in Sections IIIII and \W\ as 
well as in However, in the special case Qi = Q5, all An minimal scalars satisfy this equation. Therefore, our 
conclusions that we derived for the minimally coupled (massless) scalar field apply straightforwardly for these fields 
above. Further then, and as the main conclusion, perturbations of certain supergravity fields can drive the JMaRT 
geometries @ unstable due to the ergoregion instability Q but on the other hand, these same fields do not seem to 
produce a comparable (superradiant) instability for the D1-D5-P black holes. 

While our analysis does not apply in general to the remaining 5n — 3 scalars, we are tempted to discuss these 
in qualitative terms. Quite generally, we expect that the superradiant instability will not appear for these scalars, 
independent of most of the details of their wave equation. First, we observe that for very short wavelengths, we expect 
wave packets to propagate along the characteristic curves of the wave equation. For example, with the Klein-Gordon 
equation (flU)) . high-frequency wave packets travel along null geodesies. These characteristics are determined by the 
principle part of the wave equation, i.e., the second order terms. Now for a general scalar, we expect that these 
characteristics will in fact match the null geodesies of the Klein-Gordon field. The relevant term in the action will be 
precisely the kinetic term of a given scalar (f3"T)) . Hence we can note here that the interactions with the background 
RR two-form will not affect the characteristics. Next the background scalars will only modify the kinetic term of 
any given scalar excitation by multiplying the latter with a nontrivial overall factor. Hence the principle part of the 
resulting wave equation is only modified by an overall factor which leaves the characteristics unchanged. Hence we 
expect that short-wavelength wave packets of all of 5n scalars propagate along null geodesies of the six-dimensional 
geometry. 

Next we consider the analysis of [l(J EH which considered the superradiant instability in black strings and black 
branes in arbitrary dimensions - these were solutions of the vacuum Einstein equations. It was found that while an 
instability appears with four (noncompact) dimensions, no such instability appears in five and higher dimensions. 
Only in four dimensions did the background provide an effective potential which trapped bound states. However, it 
was pointed that in the high-frequency limit such a bound state would be following a stable circular orbit in the black 
hole background and so in higher dimensions, the absence of bound states can be related to the absence of stable 
circular orbits, as mentioned in section ITVl In the present case, we again are studying black strings in six dimensions 
or effectively five-dimensional black holes. Hence we should not expect to find stable orbits in these backgrounds. 
While a complete proof would require a new detailed analysis, the absence of a trapping potential for Klein-Gordon 
scalars certainly suggests the absence of any such orbits. The absence of such orbits can then be used to argue the 
absence of bound states and hence the absence of a superradiant instability for a general supergravity scalar. 

In passing we note that the existence of negative-energy geodesies trapped in the ergoregion of the JMaRT solutions 
can be argued on general grounds, as discussed in ||. Hence a similar reasoning to that above suggests that any of 
the supergravity scalars can initiate the ergoregion instability in these backgrounds. That is, these bound geodesies 
would correspond to trapped states in the context of a field theory analysis. The key question then becomes whether 
the corresponding scalar field modes of the field fit inside the ergoregion or whether they leak out to infinity, i.e., 
whether they correspond to a true negative-energy bound state or to a mode producing the ergoregion instability. 
The detailed analysis of § showed that both kinds of modes existed for a Klein-Gordon scalar but in particular the 
spinning modes were generically associated with the ergoregion instability. While we have not extended this detailed 
analysis to the complete collection of six-dimensional supergravity scalars, we expect that the similar results would 
be found. That is, the ergoregion instability will generically be initiated by such modes that are 'trapped' by their 
angular momentum. 

Of course, the primary motivation for the present study were the possible implications for Mathur's fuzzball proposal 
[l[ . We have already presented an extensive discussion on this topic in [H and will only comment on a few of the salient 
points here. According to Mathur's proposal, the individual microstates of a black hole are described by horizon-free 
solitons and the black hole geometry only appears after 'averaging' over these microstate solutions. Much of the 
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evidence for these ideas comes from studying certain supersymmetric solutions in five and four Q dimensions. 

However, if the fuzzball proposal is to have any substance, it must also extend to non-supersymmetric black holes. 
The JMaRT solitons provide the first family of smooth horizon-free geometries which are non-supersymmetric and 
so correspond to non-BPS microstates @. However, as mentioned in the introduction, these solutions also present 
an apparent contradiction with the fuzzball proposal. That is, the JMaRT solitons suffer from a classical instability, 
namely an ergoregion instability, and further it can be argued that this instability should be robust feature of any 
smooth horizon-free geometry carrying angular momentum [8|. However, the nonextremal rotating D1-D5-P black 
holes exhibit no comparable instability. In particular, we have shown here that these black holes do not exhibit a 
superradiant instability. Hence there is a possible contradiction for the fuzzball proposal since one would expect that 
if the ergoregion instability is common to all of the rotating non-BPS microstate geometries, then this instability 
should be reflected in the black hole geometry which is supposed to arise from averaging over the microstate solutions. 

However, this reasoning is not definitive and this puzzle still has physically sensible resolution, at least in principle. 
In particular, there are two observations which we must make about the JMaRT solitons. First, the mass and spin of 
the JMaRT and black hole solutions are in very different regimes, as described in ((HJ: M < (ai — a-2) 2 for the JMaRT 
branch and M > (ai + a2) 2 for the black hole branch. Hence the JMaRT solutions should be expected to represent 
at best a very small contribution to the microstate ensemble underlying a nonextremal D1-D5-P black hole. Second, 
the JMaRT solutions are very symmetric spacetimes. In particular, they have all of the same Killing symmetries 
as the D1-D5-P black holes, since these are simply two branches of a common family of supergravity solutions. In 
contrast, generic microstate geometries are expected to contain complex throats which do not respect these Killing 
symmetries [l[ . Hence the physical characteristics of the JMaRT solutions are likely not representative of those for 
a typical microstate in the black hole ensemble. While the typical microstate geometries should still suffer from an 
ergoregion instability, one might expect that the instability timescale becomes very long Q, especially in the 'classical 
limit' where the string coupling is taken to zero [30j. In particular, the complex throat at the core of the typical 
microstate geometries should emulate the absorptive behaviour of the black horizon in this limit, making it difficult 
to distinguish physics in these backgrounds from that in a black hole background, except on very large time scales. 
It is reasonable then that the timescale of the ergoregion instability should be a scale that grows parametrically as 
9s -> 0. 

Given that the ergoregion instability should be a generic feature of nonsupersymmetric microstate geometries, it 
would be interesting to study the dual non-BPS microstates at weak coupling for evidence of such an instability. 
As a particularly simple set of microstates have been identified to correspond to the JMaRT solitons [6(, these may 
provide a good framework to initiate such a line of investigation. For more general microstates, e.g., those expected to 
describe a near-extremal spinning D1-D5-P black hole 31], one must be careful to distinguish the expected Hawking 
radiation from radiation related to the ergoregion instability. In this regime, the latter is likely to be related to the 
'nonthermal' radiation that is expected to produce the spin-down of the black hole 32]. It may be, however, that 
when considering non-BPS configurations that the ergoregion instability provides a signature by which microstate 
geometries are more easily distinguished from their black hole counter-parts. Hence this seems a promising direction 
of research. 

Of course, another challenging problem which remains is the construction of a more or less complete family of 
microstate geometries beyond the BPS sector. While the existence of the JMaRT solitons indicates that at least 
certain non-BPS states can be described by horizon-free geometries, it is not at all clear that this property should be 
shared by all such states. However, this is certainly another intriguing research direction. 
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APPENDIX A: SIMPLIFICATION OF THE WAVE EQUATION FOR THE JMaRT SOLITON 



In this appendix we clarify the connection between the present perturbation analysis on the black hole branch of 
the D1-D5-P system and the stability analysis of the JMaRT solitons done in [8[. 

In either case, we require that the function g(r), given in eq. (0, has real roots. The JMaRT solitons then appear 
in the low-mass regime of (|SJ), M 2 < (ai — e^) 2 where < 0. While the metric may appear singular at r 2 = r+, one 
can impose a series of constraints that ensure that the solutions are free of singularities, horizons and closed time-like 
curves. This task leads to the construction of the JMaRT solitons @ • These solitonic solutions have an appropriate 
circle that shrinks to zero at the origin and the constraints ensure that this happens smoothly. First, M and R are 
re-expressed in terms of the remaining parameters - see Eqs. (3.15) and (3.20) of @, 

_ MsiCiS 5 C 5y /SiCiS 5 C 5 S p Cp 

V°^(44 c l - 444) 

2,2 C l C 5 C p + 444p / a 1 \ 

M = a\ + — ai<22 - . (AI) 

siasscsSpCp 

Then ensuring the geometry remains smooth requires imposing to 'quantization' conditions 

^ R = n 



SpCp 

a,iC\C^Cp - aisiS5S f 



R = m (A2) 



where m,n are both integers [g]. These two constraints can be put in a more elegant form by introducing the 
dimcnsionless quantities, 

/ \ 1/2 / \ 1/2 

<1, a=(™*) , (A3) 

\ a i/ \ClC 5 C p J 



with which the constraints (|Al[) can be re-expressed as 

^ ^ _. = m — n , - — ^——r = m + n . (A4) 
s + s s — s 

Again, without loss of generality, we have assumed a\ > ai > 0, which further implies m > n > 0. We also note here 
that the special case m = n + 1 corresponds to supersymmetric solutions. In this case one also has: M = 0, s = 1, 
j = 1, oi = a 2 . 

Imposing the constraints (|A2j) leaves a five-parameter family of smooth solitonic solutions. We can think of the 



independent parameters as the Dl-brane and D5-brane charges, QijQs; the (asymptotic) radius of the y-circle, R 
and the two integers, m and n, which fix the remaining physical parameters as 



isyr 
1 



~. QiQz T Q1Q5 T Q1Q5 , Kt .s 

Q P = nm r2 , = -m — , J$ = n — . (A5) 

Of course, depending on the specific application, it may be more appropriate and/or simpler to describe the solutions 
using a different set of quantities. To conclude our discussion of the JMaRT case, we note that the roots ((7j) of g rr 
can be rewritten as 

2 SiS^Sp 2 CiC^Cp 

r + = —a\ai , r_ = — a\a2 , (Ao) 

CiC 5 C p SiS 5 S p 

The wave equation in the background of the JMaRT solitons is still given by (Tl7|) . but we can simplify it by using 
the JMaRT constraints (|A1|) - (| A6|) . The results of this Appendix will make use of these constraints and so they will 
be valid only for the horizon-free JMaRT solutions; they no longer apply to the general case and in particular to the 
black hole case. Define 

y-j2 „2 „2 „2 r,2 r,2 



S1C1S5C5 

S1C1S5C5 



(A7) 
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and note that r\ — r 2 _ = aiCl2P , If and only if the JMaRT constraints are imposed, we can verify the following 
identities: 

UJ 2 M 2 UJ 2 

- -^pj -jp [~c\clc 2 p r 2 - 2s 1 s 5 s p c 1 c 5 c p a 1 a 2 + s\sls 2 p (r 2 + a 2 + a 2 - M)] = (r 2 + - r 2 _)p 2 (r 2 - r 2 + ) , 



2uj\ M 2 r r 222, 2 2/ 2, 2,2 . x\1 , / 2 , 2\ 1 2ul\ 2 2 \ n/ 2 2 \ 

-jp [~s p c p [-c^r + s 1 s 5 (r +a 1 +a 2 - M)\ + (c p + s p )sis 5 cic 5 aia 2 \ = ^-y (r + - r_)p§(r - r + ) , 



2ti)m<A Mr o , 2 2 i 

g(r) ~R l-~ ClC5C P a2 v + a i) + s i s 5S p ai(r +a 1 —M)\ 



2ujm 
9(r) 



*{r 2 + -r 2 _)pn{r 2 -r 2 + ), 



(A8) 



~g^)' ~R [ ClC5C P ai ( r2 + a a) ~ siS5S P a2(r 2 + a 2 - M)] = ( r + ~ rt)pm(r 2 - r 2 + ) , 

X 2 M 2 

- -^rr -^j [ _c i c 5 s p r2 ~~ 2sis 5 SpCiC5Cpaia 2 + s 2 s 2 c 2 (r 2 + a 2 + - M)] 

= --^(4- r 2 ) [ (r 2_ r 2 ) _ 1? 2 (r 2_ r 2 )] j 

^ f- [ Cl c 5Sp a 2 (r 2 + a 2 ) - ^^(r 2 + a\ - M)] = 2 -^±{r\- r 2 _) [-m(r 2 - r 2 _) - m?(r 2 r 2 + )} , 
^ ^ [ciC 5Sp ai(r 2 + a 2 ) - Sl s,c p a 2 (r 2 + a 2 M)] = ^± (r\ - r 2 ) [n(r 2 - r 2 ) + md(r 2 - r 2 + )\ , 

*M*L Mai a 2 = ^» (r 2 - r 2 fnm , 
g{r) ~ g{r) 

~ ^| [(? + «2)(«2 - «?) - Ma 2 2 ] = H " rl) [n 2 (r 2 - r 2 ) - ™ 2 (r 2 - r 2 + )] . 

So, in these equalities, the left-hand side was taken from (fT7|) and the right-hand side is valid only if the JMaRT 
constraints (|A1[) - (|A2[) are imposed. Hence in the JMaRT case, we leave the first line of p7|) unchanged but we can 
simplify all the terms proportional to g(r). Inserting (|A8j) into eq. (fl7|) yields the simplified version of the wave 
equation for the JMaRT solitons, 
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h = 0. 



which is exactly Equation (6.4) of Q and Equation (14) of Q. 



APPENDIX B: THE SCHRODINGER POTENTIALS 



In this appendix we give the explicit form of the function of the function 7 and of the Schrodinger potentials V± 
that are defined in (12611: 
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